Mesoscopic variations of local density of states in disordered superconductors 
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We explore correlations of inhomogeneous local density of states (LDoS) for impure supercon- 
ductors with different symmetries of the order parameter (s-wave and d-wave) and different types 
of scatterers (elastic and magnetic impurities). It turns out that the LDoS correlation function 
of superconductor always slowly decreases with distance up to the phase-breaking length and 
its long-range spatial behavior is determined only by the dimensionality, as in normal metals. On 
the other hand, the energy dependence of this correlation function is sensitive to symmetry of the 
order parameter and nature of scatterers. Only in the simplest case of s-wave superconductor with 
elastic scatterers the inhomogeneous LDoS is directly connected to the corresponding characteristics 
of normal metal. We found that in presence of pair-breaking scattering relative LDoS variations 
increase with decreasing energy. 
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I. INTRODUCTION 

Classical theory of superconductivity for impure mate- 
rials deals with average quantities. Within the BCS ap- 
proach the average fundamental characteristics of s-wave 
superconductor such as transition temperature, gap, and 
density of states are not sensitive to potential disorder. 
This statement, known as Anderson theorem, is, in fact, 
not rigid at all. In particular, it is enough to introduce in 
s-wave superconductors some amount of magnetic impu- 
rities and they suppress the transition temperature and 
gap in the quasiparticle spectrum. Moreover, in the case 
of more complex symmetry of the order parameter, even 
elastic impurities depress superconductivity. 

It is worth to note that average parameters do not com- 
pletely describe properties of impure materials, because, 
in addition, disorder induces random point-to-point vari- 
ations of all quantities. For instance, it is well known 
since 60's, that the LDoS of normal metal near an iso- 
lated impurity experiences so called Friedel oscillations 
at the atomic scaleJ* 

At the end of 70's - beginning of 80's the theory of 
weak localization was developed which described the cor- 
rections to average values of transport characteristics 
of impure electron systems caused by the quantum in- 
terference of electrons due to their multiple impurity 
scatteringf2. Even though these corrections of quantum 
origin were found to be small in comparison to the corre- 
sponding classical values, it was demonstrated that they 
have nontrivial dependences on temperature, frequency, 
and magnetic field, what makes them experimentally ac- 
cessible. Even though the quantum interference itself 
does not effect the average value of DOS, the nontrivial 
corrections to this quantity appear when, in addition, the 
interelectron interaction is taken into account. 

During mid 80's, the spatial variations of the LDoS, 
conductivity and other normal-metal properties have 
been revisited within the framework of the mesoscopic 
physicsi^ It was found that the quantum interference ef- 
fects also lead to appearance of nontrivial corrections to 



inhomogeneous characteristics, for instance, to the LDoS 
correlation function.— In contrast to the "fast" atomic- 
scale contribution of the Friedel oscillations, the latter 
manifest themselves in smooth long-range spatial behav- 
ior of the LDoS correlation function as the "slow" power 
(or logarithmic in 2D case) decay on the distances be- 
yond the mean-free path I and up to the phase-breaking 
length Irj, 3> I (each of them is much larger than inter- 
atomic distances). One can recognize the physical origin 
of such phenomenon in spirit of the qualitative explana- 
tion of the weak localization corrections given in terms 
of self-intersecting trajectories.^ 

The electron motion in impure metal has the diffusive 
character. For every pair of remote points r and r' with 
finite probability one can find the self-intersecting quasi- 
classical trajectory which starts from the point r, passes 
close to the point r', and returns back to the initial point. 
An important property is the existence of the opposite 
returning trajectory, which outcomes from the point r', 
passes close to the point r, and returns to the point r' 
following roughly the same route [see Fig. [1] (a)]. These 
two trajectories have two long joint pieces where the elec- 
tron motion is accompanied by the multiple scattering on 
the same impurities. When time-inversion symmetry is 
present, particles can move along these trajectories both 
in the same or in the opposite directions. Looking at 
Fig. [Tfa) one can see that for each trajectory there are 
the entry and the exit points of the joint routes (marked 
by circles) separated by the distances i?i and i?2 from the 
trajectory origin. Existence of such diffusive trajectories 
leads to long-range correlation of different properties, in 
particular, the local density of states. 

Quantitatively, this phenomenon can be described by 
the standard Green's function diagrammatic technique. 
The diagram describing the long-range correlations is 
shown in Fig. [Ifb). It contains either two diffusons or 
two coopcrons.^ The two-diffuson diagram describes the 
process in which the particles move in the same direc- 
tion within the two joint routes, while the two-cooperon 
diagram corresponds to the motion of particles in the 
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opposite directions. The Cooperon (diffuson) as the el- 
ement of the diagram describes the process of coherent 
scattering of electrons moving along the joint routes. The 
blocks of three Green's functions (two retarded and one 
advanced, or vice versa) are known as the Hikami boxes. 
They describe the incoherent motion of electrons in the 
domains close to the entry and exit points r and r', where 
their routes divaricate. 

Recently, STM measurements of the LDoS spa- 
tial variations have emerged as a new powerful tool 
to characterize intrinsic inhomogeneities in impure 
superconductors;^ These measurements revealed both 
rapid oscillations with typical wave vectors connecting 
characteristic points at the Fermi surface (the quasi- 
particle scattering interference patterns) and smooth 
LDoS variations. In particular, studying the oscillat- 
ing contribution provides one of the ways to establish 
the symmetry of the superconducting order parameter. 
The available theoretical description of these experiments 
is mostly based on the single-impurity approximation^, 
which becomes insufhcient at noticeable impurity con- 
centration. 

Our purpose in this article is to understand the behav- 
ior of inhomogeneous LDoS of impure superconductors 
with different order parameter symmetries at the length 
scales beyond the mean free path. We develop a theory 
which properly accounts for the collective effects appear- 
ing during coherent quasiparticle scattering on impuri- 
ties. We demonstrate that the energy dependence of the 
long-range correlation function of superconductor is in- 
deed sensitive to symmetry of the order parameter and 
nature of scatterers. The inhomogeneous LDoS of su- 
perconductor can be directly mapped on that one of a 
normal metal only in the simplest case of s-wave super- 
conductor with elastic scatterers. Presence of magnetic 
impurities or more nontrivial symmetry of the order pa- 
rameter results in the considerable complication of the 
LDoS correlation function energy dependence while the 
spatial variations in all cases do not change. 

The article is organized in the following way. In Sec- 
tion II we start our discussion introducing the Green's 
function formalism for study of the inhomogeneous LDoS 
correlation function and refresh to a reader its proper- 
ties in normal metal. Section III is devoted to study 
of the LDoS correlation function in s-wave superconduc- 
tors and it consists of two subsections treating cases with 
only elastic and both elastic and magnetic impurities. In 
Section IV we consider the problem in the case of super- 
conductor with d-wave symmetry of the order parameter 
and elastic impurities. The cumbersome technical details 
of calculation of the elements of diagrams for correlation 
function, such as Hikami blocks, cooperons, and traces of 
large number of Pauli matrices make the Appendices. 
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FIG. 1. (a) Returning quasiparticle trajectories leading to the 
long-range correlations between the points r and r'. The tra- 
jectories are characterized by two pieces of long joint routes 
within which the diffusing particles move either in the same or 
in the opposite directions. (b)Two-diffuson diagram describ- 
ing long-range part of the LDoS correlation function. Upper 
(lower) loop represents G^{E, r', r') {G^{E, r, r)) correspond- 
ingly. Shaded box represents diffuson. The two-cooperon con- 
tribution can be obtained by reversing arrows direction in one 
of the loops. 



II. INHOMOGENEOUS LDOS IN NORMAL 
METALS 

The case of normal metal represents a natural starting 
point and a convenient reference. The effects of multiple 
scattering and quantum coherence on spatial correlations 
of LDoS in disordered normal metals were considered in 
Ref. 0. It was found that two different contributions 
can be distinguished in the LDoS correlation function. 
The short-range contribution represents the Friedel os- 
cillations modified by multiple impurity scattering. The 
second long-range contribution appears due coherent dif- 
fusive propagation of normal quasiparticles. In the fol- 
lowing, we will focus namely on this long-range term and. 
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for completeness, we start with reproduction of its calcu- 
lation. 

Our purpose is to evaluate the LDoS correlation func- 
tion at the same energy 

£(-)(r-r',i?) - (iV(")(i?,r)iV(")(i?,r'))-(Af("H£;))'. 

(1) 

As the LDoS per spin is related to the retarded Green's 
function G^{E, r,r) by the standard relation, 



where I = vpT is the mean-free path. In real space 



-lm[G''iE,r,r)]/n, 



(2) 



this quantity can be expressed via retarded and advanced 
Green's functions: 

£(")(r - r', E)^^ {Re [(G«(i?, r, r)G^{E, r', r')) 



(G«(ii;,0,0))(G^(i?,0,0))] 



(3) 



-Re 



{G^iE, r, r) G"{E, r', r')) - {G^{E, 0, 0))'] } , 



where (...) implies averaging over impurities distribu- 
tion. One can use this presentation for impurity aver- 
aging within the standard Green's function approach. 
We will assume weak impurity scattering (Born limit). 
The contributions to the LDoS correlation function can 
be represented as diagrams consisting of two loops con- 
nected by impurity lines. 

The long-range contribution to the LDoS correlation 
function is given by the sum of the two-diffuson and two- 
cooperon diagrams, £^"^ = /^[jl?) +'^(2C)' Fig. [TJb). 
Both of them give the same contributions, which can be 
approximately written as 



c[f^^{r~r',E)=c[f^^{r-r',E) 
where Cn(r — r') is the cooperon. 



\B{EW 



(4) 



B^") (i;) = J dRiJ dR2G^(Ri)G^(Ri-R2)G^(R2) 



(5) 



is the "Hikami box" , and 



d^p exp(ipR) 



are the averaged Green's functions with D being the 
space dimensionality. Here r is the elastic scattering 
time, Ep — vf{p — Pf), vf and pp are the Fermi veloc- 
ity and momentum. Integration in Eq. ([5]) gives a very 
simple result 

B^") {E) = 2ttwt'^ (6) 

with V = (7Vi"^(£;,r)) as the average value of LDoS. 

Hence, the long-range behavior of the LDoS is com- 
pletely determined by the cooperon Gn(r — r') whose 
Fourier transform is well known^ 



CM 



2'KVT 



fu{lq) with fuilq) 



-r-^ for Iq <C 1, 



At r > Z 



1 / ln(/0/r) for D = 2 



TT 1 3/4r for D = 3 



(8) 



(9) 



For the 2D case the logarithmic divergency is cut off at 
q ^ l/l^, where is the phase-decoherence length. As 
the Green's functions decay at the distances of the order 
of the mean- free path, in order to evaluate the long-range 
behavior |r — r'| ^ I, we replaced in Eq. arguments 
of both cooperons with |r — r'|. 

Substituting the results ^ and dS]) into Eq. (g]), we 
obtain the long-range asymptotic expression of the LDoS 
correlation function- 

£(")(!• -r',^;) = {2uT^)^ Cl{r ^ r') 



(kpl) 



2D-2 



/2(|r-r'|/0 (10) 



with 02 = 4 and 03 = 47r^. From this result we see that 
LDoS variations are small in comparison with the average 
DoS by the parameter {kFl)^~^ . It is important to stress 
that in contrast to the short-ranged Friedel oscillations, 
they decay slowly at large distances up to l^, as \n'^{l^/r) 
for the 2D case and as for the 3D case. 



III. S-WAVE SUPERCONDUCTORS 

A. Potential impurities 

We start consideration of superconducting state with 
the simplest situation of purely potential scattering and 
s-wave symmetry of the order parameter. In this case the 
inhomogeneous LDoS is directly related to its normal- 
state counterpart. Indeed, using eigenstates expansion, 
the normal-state LDoS can be represented as 

ivi")(£;,r) = ^|^,(r)p^(£;~£;,), (11) 

i 

where ipi{r) are eigenfunctions and Ei are eigenener- 
gies of the quasiparticle states. In case of poten- 
tial scattering, the correspondin g eigenen ergies in su- 
perconducting state become ■^\/ El -I- where A is 
the gap parameter, and the two-component Bogoli- 
ubov wave function of quasiparticle state in supercon- 
ductor (C^i,±(r), l^i.±(r)) is proportional to the normal 
state wave function ([/i,±(r), Vi^±(r)) — {ui,±,Vi ,±)ipi{r), 
where Ui_± and Vi^± are coordinate-independent con- 



stants Ui 



[a 



/V2) J I + aE^/^Ef + l^'^, 



(7) 



(I/V2) ^1 - ai?,//EF+A^, |«i,±|2 + |z;,,±|2 = 14 A 
quantity commonly evaluated for superconductors is the 
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density of states for excitations which in normal state 
corresponds to symmetric combination Ncx\E,r) — 
ivi"\i;,r) + M'^\-E,r). The excitation LDoS in su- 
perconducting state, NExiE,r), can be represented in 
the form of eigenstate expansion as 



(^E,-aVE^ - A^^ 

(12) 



i,a— ±1 

\E\ 



V£^2 _ A2 



(in the second line we change the variable of the S- 
function). Such a simple connection between the nor- 
mal and superconducting LDoS is one of consequences 
of the Anderson theorem and it provides the following 
relation between the normal-state and superconducting 
LDoS correlation functions 

Cil^iE,r)/[i.i^iE)r ^ CtHy/E^^,r)/i2^r, (13) 

where (E) — 2iyE/y/E^ — is the average supercon- 
ducting DoS for excitations and Ccx {E, r) « {E, r) -|- 
£(")(_£' 1-), Even though in the following we only con- 
sider {E, r), for completeness we also present a useful 
general relation for the electronic LDoS, 



N(^HE,r)^ E mAr)\'siE-a^Ef+A 



y E±^^^Ni-\a^/E^^v) (14) 

a—±l ^ 



In particular, this general relation allows immediately to 
reproduce the single-impurity result reported for the s- 
wave case in Ref. |6|. 

Even though the long-range tail of LDoS for disor- 
dered s-wave superconductors can be immediately ob- 
tained from the normal-state result, it is useful, never- 
theless, to rederive it formally, within the Green's func- 
tion approach. This will allow us to generalize such an 
approach later for less trivial situations for which the 
above argument does not work any more. 

The long-range contribution to the LDoS correlation 
function is again determined by the diagram shown in 
Fig. [Ifb), but both the Green's functions and cooper- 
ons now have the matrix structure. In order to calculate 
them, it is convenient to use Nambu formalism decom- 
posing the Green's functions over Pauli matrices and 
the cooperons over the direct product f'' (g) f'' of them 



I^R,A _ R,A J 



c: 



is) 



k k' 
^kk'^apT^jS 



(15a) 
(15b) 



We assume summation with respect to repeated indices. 



Let us note that the Cooperon 4x4 matrix Cj,j, 
has the 2x2 block structur^>52 : 



kk' 



c 







c 



fact 



(16) 



For potential scattering the averaged superconducting 
Green's functions are given by 



,R.A 



(£;,p) 



(a«-^)2 (£;2 - A2) - £2 



(17) 



a«-^ = 1 T 



2r/B2 



and for real A they do not contain f2 component. 
In this formalism it is convenient to deal with LDoS 

is) 

for excitations TVcx which is related to the trace of the 
Green's function 



iVi^)(£;,r) = — Im[TrG^(i?,r,r) 



(18) 



Gomparing this equation with the previous normal state 
LDoS definition ^ one can see that Ni^{E,T) has an 
additional factor two since it contains both electron and 
hole contributions. In particular, the average LDoS for 
excitations is given by I'ex • The corresponding expression 
for the two-cooperon diagram can be represented using 
the Pauli-matrices decomposition: 

4 



C 



cx(2C) 



■i' ^ kk' ^ mm' ^ 



1 



Ukm = 77 Tr ( 



T T T T T-^ , 



(19a) 
(19b) 



^J^'^^J rfR25f (Ri)g^(Ri-R2).9/(R2). (I9c) 

The computational details of superconducting 
cooperon components are presented in Appendix lA II It 
turns out that its singular components are related to the 
normal-state cooperon ([7]) as 



W - E^^aA E 



C„(q) 



for i, j = 0, 1 



C^3^(q) = a(q)/2 



(20a) 
(20b) 



One can see that in the limit of normal metal, A — 0, 
the only nonzero components remained are Cqq and C^^^ . 

Computing the trace of five Pauli matrices in Eq. 
(|19bl) , taking into account (i) symmetry of B^^j with re- 
spect to indices i and j and (ii) the absence of f2 compo- 
nents in all decompositions (see Appendix [C]) , we obtain 



Ukm. = SmaBff}^ + SkoB. 



(Smk — 2(5„iO(5feo) B^Q. 



+ 2{(l~4o)4:i, + (l~<5,„o)S, 

-(1 - ^no)'Bo'2„ [Skm - 2(5feo'^mo)| 



is) 
Qkm 
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where summation with respect to the index i = 0,1,3 
is assumed in the first three terms. The components of 
the Hikami boxes are computed in Appendix lA 21 For 

i,j,m = 0, 1 



and 



2t) 



2 {E^-A^f^^{VE^-A^ + ^) 



033 



303 — 



imyEr'^ 



2t 



2r 



Due to the specifics of the cooperon clj\q) structure 
for i,j = 0,1, seen from Eq. (|20ap . we need only the 
combination Uqo + Un (A/Ef - 2 (A/E) Uqi for which 
we find a very simple relation 



Uoo + Uii (A/Ef - 2 (A/E) Uoi ^ U33 1 - {A/Ef 



and for C/33 we derive 



?733 = — 2i7rz/T 



Collecting terms, we finally obtain for the total correla- 



tion function, C, 



is) 



C 



(b) 



>c(2C) 



c 



(2D)' 



£g(r-r',£;) = -5|C/33|^ 



(§) ) [C^io^(r-r')]^ + [Cr3^(r-r')]^ 




(21) 



This result explicitly confirms the relation ([T^ based on 
Anderson-theorem arguments. 



B. Magnetic impurities 



approac h -'^^'^2 have demonstrated that the low-energy 
quasiparticle states are in fact generated for all magnetic- 
impurities concentrations meaning that, strictly speak- 
ing, the hard gap is eliminated by any amount of mag- 
netic impurities. For small concentrations, however, the 
average DoS at low energies have exponentially small tail. 

Formally, in the presence of magnetic scattering the 
Green's function can be still presented in the form sim- 
ilar to Eq. (|17p . but within the mean- field approach the 
renormalized energy and gap now are determined from 
the self-consistent transcendental equationSi^ This result 
is usually obtained taking into account the spin structure 
of the superconductive Green's functioniS, which means 
the additional increase of the matrix dimensionality to 
4x4. In the calculation of the two-cooperon diagrams 
in Nambu representation (jl9al) the number of Pauli ma- 
trices in traces is already large. That is why, for the 
sake of simplicity, we will consider below magnetic scat- 
terers as Ising impurities oriented in the same direction. 
This allows us to preserve 2x2 matrix structure of the 
Green's function, which significantly simplifies calcula- 
tions. The physical picture remains practically identical 
to the case of isotropic magnetic impurities. Even for 
such a minimum model describing the pair-breaking scat- 
tering, the calculations and results become rather cum- 
bersome. Technically similar work has been done in Ref. 
[T3I where suppression of the transition temperature by 
magnetic impurities in combination with Coulomb effects 
has been investigated. This thermodynamic problem re- 
quired summation of ladder diagrams with elastic and 
magnetic impurity lines in the Matsubara-frequency pre- 
sentation. In principle, the cooperon components at real 
energies could be obtained from the results of this work 
via analytic continuation. However, this procedure is not 
trivial at all. 

The long-range LDoS correlation function is still de- 
fined by Eq. (|19ap but the Green's functions, cooperons, 
and Hikami boxes are considerably modified by magnetic 
scattering. Detailed derivations of these objects are pre- 
sented in Appendices IB II and IB 21 Instead of Eq. (|17l) , 
Green's functions are given by 



G^'^ (p) = 



Ei-Al-el 



(22) 



In this section we evaluate the LDoS correlation func- 
tion for s-wave superconductor with magnetic scatter- 
ers. Since the seminal paper of Abrikosov and Gor'kov^ 
it is established that the magnetic impurities dramati- 
cally suppress superconductivity and strongly influence 
the quasiparticle density of states. The original mean- 
field approach (AG theory) suggested that the hard gap 
in the average DoS is not eliminated by magnetic im- 
purities. This hard gap decreases with increasing the 
magnetic-impurities concentration and vanishes at cer- 
tain critical concentration so that the gapless supercon- 
ducting state exists within small range of concentrations. 
More accurate later treatments beyond the mean-field 



with the renormalized energy and gap 



E±= E± 



A+ = A± 



i 



2n./??i-l 



(23a) 
(23b) 



where 1/tq = 1/t -|- (— 1)"/Tm, t and Tm are the po- 
tential and magnetic scattering times. Here and below 
the subscript ("— ") corresponds to the retarded (ad- 
vanced) components. The parameter ri± = E±/A± has 
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to be determined from the equation 



TmAJi^l. - 1 



E 
A' 



(24) 



The average DoS is connected with 77+ by relation^ 



(25) 



The 2x2 blocks of the cooperon Eq. ([16)) computed in 
Appendix IB II are 



C 



(sm) _ 1/47W 



1 I Xl7q\ ^0 

Tm 2r* 



T 2r*y T* \r]l - 



(26) 



and 



C 



(sm) _ 



l/47ri/ 



1 I X27q\ 

r™ 2r* y 



with 



qp 



7q 



1_ jg_U3^^7q Im7y+ , 
T 2t* / r* — 1 



_ _ ' I _9 1 I ' I 2 



(27) 



1 1 

11 2 ^2 / ' T* ~ ''IP I ^ ~ ^ 

1 + Tqp(vi=^q) / ^ \T 



Note that the component (7^™'* stih has diffusive diver- 
gency for (7 — > 0, while in C^™^ it is cut off by magnetic 
scattering. Finally, we compute the Hikami-box compo- 
nents 



(sm) 
000 



2 ElE_ 



El-Al 



Ei-Al 



3/2 



& - A2 



B 



(sm) 
klm 



4ooV+'""d, for fc,Z,m = 0,l 



with notations C+ = A_|_/_B+ = l/?7+ and 



„(sm) _ _ ' 
"fe33 ^ 



qp 



qp 



v+Cl 



(sm) _ *7rZyTqp 



El-Al 
E^C- 



772-1 



^2 - Ai 



Kji^sm; _ _ 
'^3fc3 ^ 



for fc = 0,1. 

Taking into account all these modifications, one can 
calculate the LDoS correlation function in the presence of 
magnetic impurities. Its long-ranged part comes from the 
A-block only (as it was already mentioned, the divergence 
at g — 7> in Cs is cut off by paramagnetic impurities, see 
Ea. (|B3|) ) and the final answer can be represented as 

CiTHr - r', E) ^ {2vTlf |7^|^ ^^(r - r'), (28) 

where I/tq = l/r + l/rm, Cn(r — r') is the normal-state 
cooperon (O and 



Xi + 1 1 + t/t^ 
4 1 + XiT/r,„ 

1 + _ 27?C- 1 




(29) 



E^ 



E. 



[(l + 792)Re[C+]-,?(l + C+C-)] 



with 

Vxi + 1' E^ -§r + rj+' l + T/Tm 

Normalization of the dimensionless function TZ{E) is se- 
lected by the condition Tl{E) — ^ 1 for E ^ 00. This 
function together with normalized average DoS is plot- 
ted in Fig. [2] for several values of the pair-breaking pa- 
rameter l/{TmA). As one can see, the energy depen- 
dence of the correlation function described by 'R-{E) is 
quite different from that one of average density of states. 



The most dramatic disparity is observed at small con- 
centration of magnetic impurities, for l/(TmA) = 0.1. 
In this case the function TZ{E) monotonically increases 
with decrease of energy, while v{E) goes down and fi- 
nally vanishes at E' « 0.7 A. Here, at the edge of lo- 
cal density of states gap, the function TZ{E) reaches its 
maximum. With increase of magnetic impurities concen- 
tration, when l/(T„iA) — 0.5, the function TZ{E) still 
passes noticeably above v{E), reaches its weakly pro- 
nounced maximum and remains finite when v{E) turns 
zero. The behavior of both functions becomes almost 
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FIG. 2. Examples of the energy dependences of the normal- 
ized average DoS Vcx^\E) / {2v) , Eq. H25p . (solid lines), and 
function \TZ{E)\, Eq. (|29|) . (dashed lines), which determines 
the long-range LDoS correlations for l/(rA) — 5 and differ- 
ent values of l/(rmA). Dash-dotted lines show for reference 
the BCS density of states. 



similar only in the gapless state, when l/(TmA) = 1. 
From these plots we can conclude that, in contrast to 
elastic impurities, the ratio C'^^™\E) / {v{E))'^ increases 
with decreasing energy, i.e., relative LDoS variations be- 
come stronger at smaller energy. We also found that at 
the point where the mean-field average DoS vanishes, the 
correlation function remains finite. Note again that the 
vanishing of the average DoS is not exact result but only 
a consequence of the mean-field approximation. In re- 
ality, the mean-field gap point marks the approximate 
location of transition between delocalized quasiparticles 
and the Lifshitz-tail region corresponding to localized 
quasiparticles Calculation of the LDoS correlation 
function is beyond the mean- field approach. As C{E) in- 
creases with decreasing the diffusion constant, the found 
growth of £'^™^'(£') with decreasing energy can be inter- 
preted as indication of slowing down diffusion when en- 
ergy approaches the mobility edge. On the other hand, 
our calculation is only applicable to delocalized diffusive 
quasiparticles. One can expect that in the tail region the 
LDoS correlation function should decay exponentially at 
the localization length. Therefore measurements of the 
LDoS correlation function can be used to locate the mo- 
bility edge separating delocalized and localized states. 
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IV. D-WAVE SUPERCONDUCTORS 
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FIG. 3. Schematic Fermi surface for a cuprate d-wave super- 
conductor. f„ and g„ represent, respectively, the unity vectors 
parallel to the Fermi velocity Vf and the "gap velocity" Vg at 
the n-th gap node. The ellipses illustrate the constant-energy 
surfaces near the nodes. 

Finally, let us consider the inhomogeneous LDoS for 
superconductors with d-wave symmetry of the order pa- 
rameter, which is of special interest because of its rele- 
vance to cuprate high-temperature superconductors. Be- 
havior of LDOS in these materials was extensively stud- 
ied by STM, - Formally, the long-range LDoS correlation 
function is still determined by the general Eq. (|19ap but 
the specifics of nodal gap structure has to be taken into 
account. 

We consider a two-dimensional d-wave superconduc- 
tor with electronic spectrum e{kx, ky) defined within the 
square Brillouin zone |fca;|, |fcj,| < n/a and characterized 
by the bandwidth t = e(7r/a, 0) — e(0, 0), where a is the 
lattice constant. Figure |3] illustrates the typical Fermi 
surface of a cuprate superconductor— which, as usual, is 
determined by the equation e{kx, ky) = ^, where ^ is the 
chemical potential. The order parameter of the dx2_y2- 
wave pairing state is expressed by 

Ak = Aq [cos(fci,a) — cos(fcya)] . 

The gap nodes correspond to the points k„ — ± (fcp ± fco) 
{n = 1,2,3,4), where fco is determined by e(fco,fco) = ^J■■ 
The quasiparticle spectrum is given by £k — + 
with = ^ikx,ky) — ^ and close to the gap nodes it 

can be linearized as ~ ^ (wffcf)^ + (i^gfcg)^, where k 
is the momentum measured from the node k„ and Vg sa 
Wf Ao/t^Wf . The directions of the Fermi velocity Vf and 
"gap velocity" Vg are depicted in Fig. [31 

The weak-localization effects for d-wave superconduc- 
tors were considered in Refs. [l^ and [H for arbitrarily 
strong impurity scattering. Here we limit ourselves with 
the simplest situation of weak isotropic scattering by non- 
magnetic elastic impurities. The one-particle Green's 
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function has the same matrix structure as general Eq. 



G^'^ (E) 



(30) 



Here e is the effective energy renormalized by scattering, 
7 = 7 (e) is the impurity-induced relaxation rate. Both of 
them are self-consistently determined by the self-energy 
part e ± 17 (e) = E — (e), with 



e ± Z7 (e) 



2n 2tt {i±i-f{i)) 



-2 ' 



(31) 



Nn being the number of nodes (=4 in our case), and inte- 
gration is limited by the region near one node. Perform- 
ing integration in Eq. (I3ip with the above quasiparticle 
spectrum and separating the imaginary part of T,^ (e) 
one finds the transcendental equation for determination 
of the relaxation rate 7(e): 



In- 



An 



e e 2'KViVc, 

— arctan — — 



^l^ + 72 ' 7 7 "^r. 

At zero energy it reproduces the known resulti^"— , 

2'KViVc, 



(32) 



7(0) = 70 = Ao exp 
while at large energies 

7(e>7o) ~ -^1 



2 ln(g/7o) + 1 
The real part of S^j" determines the value of e: 




123456789 10 

E [in units Yo/ln(A(,/yo)] 



123456789 10 

E [in units Yo/inlAo/yJ] 



FIG. 4. (a)The energy dependences of the parameter e and 
the relaxation rate 7. (b) The energy dependence of the func- 
tion Fd which determines the LDoS correlation function in 
Eq. (EHI 



^fefc^^ in the Pauli- matrix expansion, Eq. (jl5bp . are sin- 
gular and these singular components are connected by 



relation C, 



(sd) 
00 



^11 



-c. 



(sd) 



(sd) 
33 



Therefore 



the matrix structure of d-wave cooperon is different from 
structure of the s-wave cooperon, Eq. ([20]) . Derivation 

of the singular cooperon CAisi) = C'oo'^ ''('?) for arbitrary 
energy is presented in Appendix ID II and result can be 
presented as 



Cd(g) 



7 



(35) 



where parameters 7, Ncx\ and the diffusion coefficient 
V are energy dependent with 



2'KViVf, 



fin- 



Aq 



7 arctan - 



E. (33) 



V{E) 



1 



- arctan - 

7 7 



27 In 



- arctan - 

7 7 



(36) 



For small i?, e = ln(Ao/7o)i?. When the energy is large 
{E » 70) 



E 



ln(Ao/7o) 
ln(e/7o) 



Note that our result for Cd{q) is smaller by factor four 
than the result of Refs. [H and[l^. The origin of this dis- 
crepancy is discussed in Appendix ID II For more trans- 
parent representation of the energy dependence in Eq. 
([55]) . we note a useful relation 



As follows from the structure of Eqs. ([5^ and the 
energy dependences of e and 7 have scaling form, e = 
7ot^£[ln(Ao/7o)i?/7o], 7 = 7o^?7[ln(Ao/7o)£'/7o]- These 
dependences are presented in Fig. HJ^a) 

The LDoS in superconducting state is determined by 
the integral of the imaginary part of the Green's function 
([50)1 which givesii 



7 In 3 



^0 



g2 ^ ^2 



2e arctan - ) . (34) 
7/ 



The structure of cooperon for d-wave superconductors 
has been investigated in Ref. [l^ and It was demon- 
strated that only diagonal components of the cooperon 



Ni^HEME) = N^^\0)Vmi + (6/7) arctan(e77)]. 



Calculation of the Hikami boxes b'i'^j defined by Eq. 
(|19c|) with Green's functions ([50)) (see Appendix ID 2p . 
results in: 



7(d) 
-"ooo 



"101 



VfVg 



- «7]^ [e - 



(d) 
303 



n r~ 

e - 



-Z7]J(°), 
Z7] J(2), 



VfVg 

"330 ^ '-'033 ^ '-'oil ^ "110 ^ 



(d) _ R(d) _ R(d) _ R(d) _ Nn_ ^(2) 
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where 

J(0) 
J(2) 



IGtt (e + 17)^ £7 



(e + ^7) e 
^ ^ arctan — 

£7 7 



327re7 



(e — 17) e 

1 H arctan — 

£7 7 



Further summation of these, non-zero, components of 

^inj ^1- CSEI with corresponding traces of five Pauli- 
matrices products as the coefHcients allows to present the 

tensor Ukm hi the diagonal form uj^!^ = Uk^km with 



C/0;2 



(d) 
000 

(d) 
000 



2B 



(d) 
101- 



In result, the LDoS correlation function Eq. (|19al) for the 
d-wave case takes the form 



4d,)(r-r',£;) = 4 



(d) 
000" 



^"101 



7(?) 



C2(r-r',ii; = 0), 



(37) 



with 7 (e) and e{E) have to be determined from Eqs. 
and (iMt. 



Fd(x) 



1 + 1/a 



(1 + X arctan x) I 2 



1 



■ arctan x 



'2 arctan'^ x} 



[2 for X = 



\ 3/2 for X > 1 ' 
and the zero-energy cooperon in real space is given by 

Cd(.,i^ = 0)^|^^ln^. 

Using this result, we can rewrite LDoS correlation func- 
tion in a more transparent form 



327r2 (t;2) 



■In" 



r — r' 



(38) 



We can see that, in contrast to s-wave superconductors 
with potential impurities, the energy dependence of the 
correlation function is not determined by the average 
density of states. Additional dependence characterized 
by the function Fd[e{E)/'y{E)] appears. Both the ra- 
tio e/7 and function Fd have universal dependences on 
the scaled energy (i?/7o) ln(Ao/7o). The energy depen- 
dence of Fd is plotted in Fig. lUJb). We observe the same 
tendency as for an s-wave superconductor with magnetic 
impurities: the relative LDoS variations become stronger 
at smaller energies. 



V. DISCUSSION 

It is known that observable quantities of disordered 
metallic systems exhibit large mesoscopic fluctuations 1^ 
The large-scale correlations of the LDoS demonstrates 
long tails extending up to r ~ Z^;^- As we have seen 
above, due to diffusive propagation of quasiparticles, such 
property of the LDoS correlation function remains valid 
also in superconductive state. This correlation function 
for all cases depends on scatterers potential V in the 
same way: C ~ and its spatial dependence is deter- 
mined only by the dimensionality of superconductor. Our 
analysis demonstrates, however, that the energy depen- 
dence of the large-scale LDoS correlations carries valu- 
able information about the order parameter symmetry 
and character of scattering in various superconductive 
systems. Indeed, as we have seen above, while the spa- 
tial dependence of £{r — r',E) is the same for all types 
of considered superconducting systems and is determined 
by the square of normal-metal cooperon, the energy de- 
pendence of the magnitude of such correlator can serve as 
the fingerprint of superconductor intrinsic properties. In 
the reference case of s-wave superconductor with elastic 
impurities such energy dependence is just given by the 
square of the BCS quasi-particle density of states. More 
complex symmetry of the order parameter manifests itself 
in significant changes of the energy dependence of LDoS 
correlation function. Investigated above the case of d- 
wave pairing (see Fig. UJ^b)), demonstrated that the cor- 
responding LDoS correlation function, normalized on the 
appropriate square of the LDoS, already depends on the 
quasi-particle energy, although this dependence is rather 
smooth. 

To illustrate the role of pair-breaking scattering, we 
analyzed the case of s-wave superconductor containing 
magnetic impurities. It is worth to discuss our results 
in light of the extended AG theory in Ref. [ll|, which in- 
vestigated the low-energy behavior of the average DoS 
in the framework of the nonlinear sigma model. The 
authors of Ref. O demonstrate, that AG results, estab- 
lishing the formation of the hard gap in the local density 
of states, correspond to the saddle-point solution of the 
proposed effective action. Their non-perturbative exten- 
sion beyond the mean-field approximation results in ap- 
pearance of sub-gap exponential tails in the density of 
states, what should lead also to appearance of nonzero 
moments of this physical value for energies below the AG 
gap. Plausibly, our perturbative calculus of the second 
moment of local densities of states C^^'^^r — r',E) (see 
Fig. [2|) indicates the same phenomenon. Our analysis 
allows for a straightforward generalization to other simi- 
lar situations, such as interband scattering in multiband 
superconductors with different signs of order parameter 
in different bands {s± state), a widely discussed model 
for iron pnictides and selenides. 

Modern STM technique, in principle, allows to probe 
the long-range correlations studied theoretically in this 
article. The challenge is to perform scan over large ar- 
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eas with sizes exceeding mean-free path. Another inter- 
esting opportunity provided by STM measurements of 
the LDoS correlation function is to extract and study 
the temperature-depending phase-breaking length l,p{T) 
of quasiparticle excitations in superconductors of differ- 
ent nature. As far as we know, such a problem was never 
considered neither theoretically nor experimentally. 
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Appendix A: s-wave superconductor with elastic 
impurities 



To proceed, we represent the matrix impurity line as 



i,j=0 



(A2) 



with li,- 



^ p(33) 



where, following Refs. I2llandl22l we introduced notation 
for the trace of four Pauli matrices 

We will need only the following components of E^/ 
which have a simple block structure2i*2^ 



Em 



fO 0' 

fO 

fi 0' 

f2 



^(11) 

^(10) 







;.0 





-f2 



f33) 

Since E^j is diagonal, we have Uij — ^aiSij with ao = 
ck3 — ~ai — — a2 — 1. 

We use a similar decomposition for the polarization 
operator 



1. Calculation of s-wave cooperon 





P Pi 



FIG. 5. Graphic presentation of Eq. (|A1|I for the supercon- 
ducting cooperon. 

In this Appendix we consider the cooperon structure 
for a s-wave superconductor with weak potential scat- 
tering. The superconducting cooperon has four Nambu 
indices and obeys the equation graphically represented 
Fig. El 



TTTTl 



-i7icti 



(Al) 



where we again assume summation with respect to 
repeated Nambu indices. Here the superconducting 
Green's functions G^'^ averaged over impurities are 
given by Eq. (jl7p . For Born potential impurities, the 
matrix impurity line is given by 



aai _ -3 -3 
11 a a. — U-T^^. To 



with u = UiV^ — 1/{2'kvt). 



711/1 



(A3a) 



with 



-Tr 



2iy 



= eI''^p,, 

Irn ^3 



■jjUpTiLTp_qT„ 



(2^) 



(A3b) 



and 



where 



P,,(q) = ^1 dCpG^pGj^p-q^^ , (A3c) 



. .) p means averaging over the Fermi surface. 
Nonzero components of Py are only for subscripts (i, j) = 
(0,0), (1,1), (0,1), (1,0), (3,3) and they can be straight- 
forwardly calculated as 



-Poo = 7rr7q 



E' 



E^ 



Pn = TTTJq 



A2 



E^ - A2 



Poi = Pw = 7rr7q 



EA 



E^ - A2 



, -P33 = 7I'T7q 



with 

7q 



1 + t2 (vf q)' 



1 - ( (Vf q) 



(A4) 



Due to the block structure of the corresponding compo- 
nents of Ef^ , the 4x4 matrix S is composed of two 
independent 2x2 submatrices, [0, 1] and [2, 3] blocks 



S = 



Sa 0' 
Sb 



(A5) 
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These blocks can be explicitly found 

Sa = r° (Poo + Pii) + f^P33 + 2f iPoi, 
5b =fO (Poo -Pii)-r3 P33. 

Note that, due to the relation Pqo — Pn = P33, the com- 
ponent 5*22 vanishes. 

Equation for the Pauli-matrix components of the 
cooperon cfj^ = jTj^a^pp^T^p^^^ now takes the form 



This allows us to present the cooperon components in the 
following form 

C^^ ^C^{~A/Ey+', (z,.7 = 0,l) 

with 

'-'00 ~ £;2 _ ^33 ■ 

Note that in the normal-state limit, A — ^ 0, only the 
components Cqq^ and Cgg'' remain singular. 



2. Calculation of s-wave Hikami boxes 



The matrix C{s) has the same block form as S meaning 
that this 4x4 system splits into two independent 2x2 
subsystems which allows us to obtain analytical results 







2 ^2 _ ^2 



(s) _ ;^ j -1 

^ "2\0 (l-7q)-' 



(A6) 
(A7) 



We can see that C33''(q) equals half of the normal-state 



cooperon. 



C^3Hq) = oCn(q) 



D 



(A8) 



The singular part of C^"* for g — is 



C 



u/2 



(vfq)' 

u/2 



EA 



-J3 



7 To 



2£;2-A2 " E^-A' 



(vfq)' 



EA 



EA 
E^-A^ 



(A9) 



Hikami boxes B^^j defined by Eq. (|19cp . In k-space B^^j 



In this Appendix we present details of calculations of 
kami boxes B^^j defined 
are given by the integrals 



B 



where ' are defined in Eq. ()15ap . 

For example, the component /Sqoq is given by the inte- 
gral 

4t = ^ / d^-r ^—^T-r : 

{a'^f {E^-A^)-e (a^f {E^-A^)-e 

Performing integration and substituting a^'^, see Eq. 
(IT71) . we obtain 



Z3 5,^(ek)g^(ek)<?/(ek) 



(2^) 



(AlO) 



(s) __nnyT 
"000 ^ 



E^ (VP2 - A2 + 3V2r) 



2 (£;2- A2)3/2(/B2~A2 + z/2r) 



As = {A/ E)g^^'^ , the components S-^^- for i, m, j = 
0, 1 are connected with the Sooo by the simple relation 

fiii-=S^oUA/i?)^+^"+". 

Two remaining nonzero components, ^033 '^3031 
given by 



R,A 



,(s) 
033 

303 



E / 



e - (Vi?2 - A2 + ^) e - - A2 - ^ 



r 



and evaluation of the integral gives 



B\ 



033 ~ 



ittvEt'^ 
" 2V-B2 - A2 ' 



303 033 2V£^2 - A2 E"^ - A^ + ^' 



Appendix B: s-wave superconductor with magnetic 
impurities 

1. Calculation of cooperon 

In this appendix we present computation details of the 
cooperon for s-wave superconductor with Ising magnetic 
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impurities polarized along z axis. The equation for the 
cooperon graphically presented in Fig. [S] is again given 
by Eq. (jAip but the impurity line here is now composed 
by potential and magnetic contributions 



with Ura = l/(27ri^rm). In order to obtain the explicit 
expression for the cooperon we follow the same route as 
in the case of potential impurities. We again present the 
impurity line in the form \A2\ . where the 4x4 matrix 
Uij is now given by 



U^ 



uE. 



(33) 



Using this presentation and formulas for the polarization 
function (see Eqs. (|A3aP and (jA3bp ) one can write the 
equation for the cooperon: 

Cf;^^ = ^5.,+i^uAiC'i^^\ (Bl) 

To evaluate the matrix Sn^ we again have to compute 
the integrals defined by Eq. (|A3cP but accounting for 
scattering from magnetic impurities. For nonzero com- 
ponents, we obtain 



-Poo 



3 7q- 



9 



Pii 
P33 



1 



3 7q- 



3 7q 



where 77+ is introduced by Eq. (j24p and 7q can be ob- 
tained from 7q defined in Eq. (|A4[) by the replacement 
of T with the energy-dependent quasiparticle relaxation 

time Tqp = 1/(2 ImWE'^ — A^). The useful expression 



for the latter parameter is derived in Appendix IB 31 



1 



Xi 



1/Tqp = - H with xi = 



1 



1 



Equation (jBl[) for the cooperon again splits into two 
independent 2x2 blocks 



^(sm) _ 1 



where the blocks of the matrices u and S can be evalu- 
ated as ^ (u/2)f3 4- (u^/2)f°, = - (u/2) -f 
(Mm/2)f", and 

= r° (Poo 
S^^T^ (Poo 
Deriving relation 



^qp7q 



1 

27 ' 



2fi RePoi, 
2if2 ImPoi. 



Xi 

^ ' in 

Re?7-) 
I ? 



Xi 
2t 



1 

2^ 



we can formally present C^^™'' as 



(sm) 



To 



1 ''qp7q 



2t 2T7T1 



Re?7+ 

'■'■qp7q"|7:T 



1 



-T^ H 



-T3rqp7q 

, -1 



Xr_^A 
2t 2t„ 



1 



(uf ^ -I- Umf°) 



After some algebra this expression can be transformed to 
much simpler form 



C 



(sm) _ 



l/47rzy 

l-7q 



1 , 'TqpT'qXl ( \ \ 

r2 ^2 



.2 ^2 



'qp7q I 2 



2 

Re 77^ 



(B2) 



Similarly, in order to compute the block cjj™'' , we de- 



rive relation 



3 7q 



2r 2t„ 



X2 i 
2t 2t„ 



r T, 



1 



Im 



with X2 = (h+P ^ 1) / |^+ ^ l|- This allows us to 
present solution for C^™'' in the form 



(sm) 



c 



+ T-qp7q 



X2^J_ 

2t 2t„ 















f3 + 


(7- 



Im 77_|_ 



After straightforward algebra one can reduce it to 



^(sm) 



l/47r!y 



l-Tqp7q(l/T + X2/-rm) 



X2Tqp7q ( \ 1_ 

r2 ^2 




1 '^qp7q ( \ L 

-2 _2 



+«Tqp7qr 



(B3) 



2. Calculation of Hikami boxes 



Let us calculate Hikami boxes (lAlOp in the case of s- 
wave superconductor with magnetic impurities, i.e. using 
the matrix Green's function (22). For the component 



Pqq™ one finds 
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(sm) 
000 



e-i^ejE^-^X + it 



with 



-, 2 r 

— (1 + ik) — (1 — ik) 



This integral can be easily evaluated for relevant values 
of ^: 



ITT I (1 + 3'"^) 

8k2(1 +m) I 1, for A* = 2 



T+wrTTT+^' tor ^ - U 



Re^/^;2 -A^ 



In result 



(sm) 
000 



qp 



3/2 



77_ 



1 



For the convenient representation of the remaining 
components Bimj we introduce notation C+ = ^+/E+ — 
1/77+ . When the subscripts are set to k,l,m = 0, 1 one 
finds 



B 



(sm) 
kml 



El 



(El - K\ 
Other non-zero components are 



3/2 



^2 _ A2 



p,(sm) 
-°033 
(sm) 



B. 



303 



1/ / 



Rc./£;2 - A 



Rc^AP 



2Ta 



E± 



and 



p,(sm) _ „ 
Bi33 - C+i'l 



Rc^/^2 _A2^ 



;B^^\ Si^-U CIB^^^ for * = 0,1. 



^2 _ A^ 



I 

following relation for the imaginary part 

1 



Vi = 



r„,,A 



■Re 



3. Quasiparticle relaxation time Tq, 



The quasiparticle relaxation time l/rqp(£') 



^4 



1 



I r]r Re J?/2 - 1 



rji Im 



TmA W+ - 1| 

2 Im J El — A^_ determines the value of diffusion con- 
stant of quasiparticles in presence of magnetic impuri- 
ties. In this Appendix we derive a useful formula for this Using this result, the factor Im Jrjl ~ 1 in the right-hand 
parameter. From the definition of A-|_, Eq. ()23bp . one 
obtains the relation 



side of Eq. (|B4[) can be transformed as. 



l/rqp(i;) = 2Im 



A«,/772 -1 



= 2AIm Jr]l - 1 + - - — . (B4) 

Presenting the parameter ry+ as the sum of its real and 
imaginary parts, rjj^ = rjr -\- irji, we find Eq. (|24p from the 



ImA/?72 _i = . 



I ?7,+(Im 



1 



2t™A 



TmA 



1^1 -1| 
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Substituting this result into Eq. (jB4[) . we finally obtain 
presentation 



(B5) 



which was used in Appendix IB II 

Appendix C: Trace of five Pauli matrices 



First of all let us symmetrize trace (|19b[) with respect 
to two indices: 

— ir (^r r r r r-* j = — lt[[t-' t + r t-* j t r r J. 

Let us recall that 



s=l 



+ (1 - J,o) (1 - 4o) 
and write down the trace of three Pauli matrices: 

{siks and Soknm here are 3D and 4D Levi-Civita symbols). 

A further step is the calculation of the trace of four 
Pauli matrices: 



1 



^ ' r T r 



s=Q 

= (1 ~ ^jo) [SjkSnm + Sjm^kn — Sjn {Skm — ^SkoSmo) 

nm 

We used the relation 



d-joT ) T 



s=0 

and the expression for product of Levi-Civita symbols 

3 

eOj/tsEOsnm = (1 ^ — ^ko) {Sju^km — ^jm^kn) ■ 

s=0 

Finally one finds the required symmetrized trace of five 



Pauli matrices: 

^ X^k 'im x^j , ^? ^/c^^n 'im 'lil 

— ir T T T T T-' + r T T T T \ — 

4 

— Sij {SmoSkn + SkoSmn + SjioSmk — '^SmoSkoSnO + i^Qkrim) 
+ SiO (1 — Sjo) [SjkSnm + SjmSkn ~ Sjn {Skm ~ 2(5feo<5,„o) 
~r 1^0 jnm 

+ SjO (1 — ^io) SikSnm + SimSkn — ^in (Skm ~ 2(5A:0'5mo) 

Appendix D: d-wave superconductors 
1. Calculation of d-wave cooperon 

The cooperon corresponding to propagation of quasi- 
particles for superconductors with d-wave symmetry of 
the order parameters can be found from the same Eq. 
(|Aip as it was done above. We again employ its expan- 
sion in terms of Pauli matrices (see Eq. (jlSbp V The anal- 
ogous expansion IIJIJ^ = n,„„r™^T"^^ we apply also to 
the polarization function, the only nonzero components 
of which are IIoo, IIii, and H^. 

According to Refs. llSlandllaonlv the diagonal compo- 
nents of such cooperon are singular. Keeping only these 
components, we find that they obey relatively simple lin- 
ear system of equations: 
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In accordance with Refs. [la and|Ia, the ansatz 
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(sd) _ ^(sd) _ 



00 
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(sd) _ ^_(sd) ^ 



33 



Cd (D2) 



reduces the right hand side of all four equations to 
(1 — unii)Cd. As we will see below, nMa^q,^ — )■ 1 for 
g, f2 — meaning that the relations (|D2I) correspond to 
the singular eigenvector. Taking the linear combination 
of equations (iPTal) + (IDlbD + (|DTc)) - (|Dldp , we obtain 



(1 - un,,) ( c, 

which gives 



(sd) 

00 



(sd) 
33 



a 



(sd) 
22 



Cd(g) = («/4)/[l-un,,(g,0)] 



(D3) 



We noticed, however, that this result is 4 times smaller 
than the one reported in Refs. [l^ and 
presentation of derivation given in Ref. 



16 



A detailed 



16 allows us to 



trace the origin of this discrepancy. The authors of Ref. 
[161 obtained their expression for Cd by substitution of the 
relations (|D2p between the singular parts of c'ff^ into 
Eq. ([DTbj) . what led to Cd = u/ [1 - -un,i(g, 0)]. This 
step, however, is problematic since side by side with the 
singular parts, c'f^'^'' also contain regular contributions: 

^ti^^ = ^^'^''d'^^ + "^ith a2 = -1; a; = 1 for i 7^ 
2 and c; are constants. Substituting such presentation 
into Eq. (jDlbp , we immediately see that the constants ci 
contribute to the nominator of Cd. As follows from the 
result (|D3p . this contribution amounts to its four times 
reduction. 

Let us pass to evaluation of the trace of the polarization 
function: 



15 



n..(g,17) 
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e-fi/2-i7) -e 
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with VL = de/dE 57. Note that the vanishing of 1 — 
ullii{q, fi) for g, 17 — !• leading to diffusive behavior fol- 
lows from general identity 
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2i7 



(D5) 



with energy-dependent diffusion coefficient 
/.,2\ 1 -t- - arctan - 



ViE) 



■ 7 



27 In 



An 
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■ arctan ■ 



(D8) 



Its value for zero-energy was first obtained in Refs. [if 



, . 1 T, J ■ J r ii, u J c V c and 16. Substituting Eq. (ID7I) into Eq. (ID3I) and using 

which can be derived from the above denmtion of , . , ^ \7i ^ n , . 

-j-j / relation 7 = [tt / Z)uN i^^^d valid lor all energies, we obtain 

^1^, ' ., r ,1 ■ r,i the presentation of the cooperon in Eq . (1351) of the main 

Ihe next step is to perform the expansion ol the po- i^g^t ^ — 

larization operator over q up to quadratic term; 



un,,(g,0)-l = - 
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- arctan - 

7 7 



472 In -,^2 — + 1 arctan 



(D6) 



with (u^) = (wf -I- w|)/2. For the trace of the polariza- 
tion operator at finite frequency and g = 0, n^,(0, 17), we ^^^^^^ parameter: 
derive the following relation 



2. Calculation of d-wave Hikami boxes 

In this appendix we calculate the Hikami boxes for 
two-dimensional superconductor with d-wave symmetry 



Mn,,(o,ri) 



^+{E+%)-^-{E-^) 



As ReS^ = - ^, the difference T.+ {E+%)-Y.^{E-^) 
can be represented as 



^+2 



^-2 



2^7 



= f7 - 17 - 2i7, 
leading to a very simple result for small 17 

n,j(0,17) w 1 + il]/27. 
Collecting terms, we obtain 

un,,((7,17) = 1 + -L (il7 - 
27 



(D7) 



Bf2j= J dUi J dR2Gf(Ri)G;^(Ri + R2)Gf(R2) 



Gf(k)G;^(k)Gf(k). 



The Green's function for d-wave superconductor is deter- 
mined by Eq. pop . which we rewrite in the form 



with 
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Hence the block B)'^ can be rewritten as 

tin J 



^(d) ^ JV^ 
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with a = e + 17. The only non-zero components are with fi = 0,2. Then the blocks Bimj are expressed as 

g(d) g(d) „(d) „(d) „(d) g(d) „{d) , , 
''000'''011'''l01'''ll0'''033'"303'"330 '^"'^ ''^^^ 

plicit expressions for them one has to carry out the inte- 
grals: 
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The integrals J^^^ can be computed explicitly as 
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